The Vortex State in a Strongly Coupled Dilute Atomic Fermionic 

Superfluid 
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We show that in a dilute fermionic superfluid, when the 
Fermions interact with an infinite scattering length, a vortex 
state is characterized by a strong density depletion along the 
vortex core. This feature can make a direct visualization of 
vortices in Fermionic superfluids possible. 

PACS numbers: 03.75.Ss, 03.75.Kk, 03.75.Lm, 03.75.Hh 



The existence of stable vortex states is one of the most 
spectacular manifestations of superfluidity in both Bose 
and Fermi systems. In Bose dilute atomic gases isolated 
quantized vortices [1] and subsequently large arrays of 
vortices [2] have been observed already. The quest for the 
observation of superfluidity in dilute Fermi gases started 
as soon as the first experimental evidence of an atomic 
dilute Fermi degenerate gas was published [3]. In Bose 
systems methods for establishing unambigously both the 
superfluidity and the quantization of vortices exists [4], 
while most of the suggestions made so far for Fermi sys- 
tems are rather indirect. The presence of vortices in BEC 
was confirmed by analyzing the density variations of an 
expanding BEC cloud [2,4]. The Bose-Einstein condensa- 
tion (BEC) creation was confirmed by studying the char- 
acter of the expansion of the atomic cloud after the trap 
was removed [5]. A recent experimental result in Fermi 
systems [6] suggests that a hydrodynamic expansion of a 
Fermi superfluid is a plausible scenario [7] . 

In Fermi systems significant density variations due to 
the presence of vortices are not expected [8]. We have 
reasons to expect however, that under certain conditions 
the density variations induced by the presence of one or 
more vortices could be akin to those in the case of BEC, 
as a similar feature was recently established in the anal- 
ysis of the spatial structure of a vortex in low density 
nuclear matter [9] . At densities significantly smaller than 
nuclear saturation densities, the superfluid gap in homo- 
geneous neutron matter can attain values rather large by 
normal standards, A 0.25£p, where A is the value of 
the gap and Ep is the Fermi energy. Even though the 
gap is still smaller than the Fermi energy, such values 
proved sufficient in the case of low density neutron mat- 
ter to lead to major density depletions in the vortex core. 
As Pitaevskii and Stringari note [6], the observation of 
quantized vortices in a dilute Fermi gas would provide the 
ultimate proof that the system has undergone a transi- 
tion to a superfluid state. A vortex is just about the only 



phenomenon in which a true stable superflow is created 
in a neutral system. Other phenomena would only some- 
what indirectly be affected by the onset of superfluidity. 

Recently, a new ab initio calculation of the properties 
of low density nuclear matter became available [10]. A 
particular result of this analysis, which is valid in princi- 
ple for any Fermion system, concerns the properties of a 
two fermion species interacting with an infinite scatter- 
ing length. These authors have shown that the energy 
per particle of such a normal Fermi system is 
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(1) 



with a at ~ 0.54. This result was obtained for a sys- 
tem of fermions interacting with a short range attractive 
potential with a zero-energy bound state. As long as 
kpro -C 1, where ro is of the order of the radius of the 
potential, one expects on general grounds that the energy 
per particle of such a system is proportional to the energy 
per particle of a non-interacting system. The scattering 
length a and the effective range ro parameterize the low 
energy behavior of the s-wave scattering phase of two 
particles fccot<5(fc) w — 1/a + rok 2 /2, where k is the wave 
vector of the relative motion. A second result of this ab 
initio calculation concerns the energy per particle of the 
superfluid phase of such a system [11], namely that 



5 



(2) 



with as ~ 0.44. A simple estimate of the correspond- 
ing value of the superfluid gap, using for the conden- 
sation energy the weak coupling BCS value — 3A 2 /8£f 
leads to A ss 0.4e^. Such large values for the gap A are 
not compatible with the BCS weak coupling limit, when 
fcp|a| <§; 1 and a < 0. It is well known that in this limit 
the BCS approximation [12] leads to a too large value 
of the gap and that polarization corrections lead to a 
reduction of the gap [13], namely to 

A= (£) ^ exp (^ 

A recent analysis [8] of the vortex state in a dilute super- 
fluid Fermi gas, using the simple BCS value for the gap 
(which exceeds by a factor of sa 2.2 the true gap value) 
shows a relatively modest density depletion in the vortex 
core of about 10% at most. 
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The possibility that the value of the superfluid gap can 
attain large values was raised more than two decades ago 
in connection with the BCS — > BEC crossover [14,15]. 
One can imagine that one can increase the strength of 
the two-particle interaction in such a manner that at 
some point a real two-bound state forms, and in that 
case a — > — oo. By continuing to increase the strength 
of the two-particle interaction, the scattering length be- 
comes positive and starts decreasing. When a 3> r 
and a > the energy of the two-particle bound state 
is 62 w —Ti 2 /ma 2 . A dilute system of fermions, when 
nr^ <C 1, will thus undergo a transition from a weakly 
coupled BCS system, when a < and a = 0(tq), to a 
BEC system of tightly bound Fermion pairs, when a > 
and a = O(ro) again. In the weakly coupled BCS limit 

the size of the Cooper pair is given by the so called co- 

h 2 ad- 
herence length £ oc m £ , which is much larger than the 

interparticle separation k \ F = 2-k /kp. In the oppo- 
site limit, when kpa <C 1 and a > 0, and when tightly 
bound pairs/dimers of size a are formed, the dimers are 
widely separated from one another. Surprisingly, these 
dimers also repel each other with an estimated scatter- 
ing length w 0.6. ..2a [12,16] and thus the BEC phase is 
also (meta)stable. The bulk of the theoretical analysis 
in the intermediate region where k F \a\ > 1 was based 
on the BCS formalism [12,14,15,17,18] and thus is highly 
questionable. All these authors have considered only the 
simple ladder diagrams in the particle-particle channel. 
The inclusion of the additional " bosonic" degrees of free- 
dom [18], which represent nothing else but a two-atom 
bound state, falls into the same approximation scheme, 
which includes only ladder diagrams in the atom-atom 
channel. Even the simplest polarization corrections have 
not been included into this type of analysis so far. In par- 
ticular, it is well known that in the low density region, 
where a < and k F \a\ <C 1 the polarization corrections 
to the BCS equations lead to a noticeable reduction of the 
gap [13]. Only a truly ab initio calculation could really 
describe the structure of a many Fermion system with 
k F \a\ 3> 1. In the limit a = ±oo, when the two-body 
bound state has exactly zero energy, and if k F ro <C 1, 
one can expect that the energy per particle of the sys- 
tem is proportional to s F = h 2 k F /2m, as it was recently 
confirmed by the variational calculations of Refs. [10,11]. 

As in any other Fermion systems [19], the knowledge 
of the energy per particle as a function of density in the 
homogeneous phase, Eq. (1), allows us to construct the 
normal part of the Local Density Approximation (LDA) 
for the energy density functional (EDF). The extension 
of the LDA to superfluid systems presented in Refs. [20] 
in conjunction with Eq. (2) lead us to a rather well de- 
fined EDF, appropriate for Fermion systems with infinite 
scattering length, namely 



£{r)n{r) 



— r(r)+Mr) 5 / 3 + 7^ 



. (3) 



n(r)=J> Q (r)| 2 , r(r) = £ |V^(r)| 2 , 

a a 
a 

Here n(r) and v(r) arc the normal and anomalous 
densities and r(r) is the kinetic energy density, all of 
them expressed through the quasiparticle wave functions 
(u a (r) , v a (r)) . The summation over a should be in- 
terpreted as a sum or integral when appropriate, over 
all quasiparticle states with the Bogoliubov-de Gcnncs 
eigenvalues E a > 0. The dimensionless couplings m/m*, 
(3 and 7 have to be chosen so as to reproduce in the case 
of infinite homogeneous matter the expressions (1-2) for 
the energy per particle. Unfortunately the parameter- 
ization of the EDF is not unique and its form has to 
be restricted using some additional considerations. One 
possible choice corresponds to m/m* = ajv, P = and 
7 w —6.64. We shall refer to this case as the parame- 
ter set I. Another possible choice of parameters, referred 
to as the parameter set II, corresponds to m/m* = 1, 
(3 = 3(a N - 1)(3tt 2 ) 2 / 3 /5 and 7 w -7.20. Since the 
anomalous density v(r) is a diverging quantity, a regular- 
ization procedure of the above expression (3) is required. 
This was performed according to the formalism described 
in great detail by us elsewhere [20] . The renormalization 
procedure amounts to replacing 7 with a new well defined 
running coupling Jeffif)- Even though in this renormal- 
ization procedure there is an explicit energy cutoff ap- 
pearing in the formalism, no observable is affected by its 
presence, when this energy cutoff is chosen appropriately. 
Up to the overall factor h 2 /m, the energy density in Eq. 
(3) has the overall scaling £(r) oc [n(r)] 5 / 3 , which is ex- 
pected from dimensional arguments for a system with an 
infinite scattering length and kpro -C 1. In principle, 
there are an infinite number of possible local energy den- 
sities satisfying this requirement. In particular, one could 
have considered for the superfluid part of the functional 
for example a contribution of the form cx |^(r)| 5 / 3 . Since 
a term cx |^(r)| 2 works in both BCS and BEC limits 
and since the emerging pairing field A(r) has the sim- 
plest form and the renormalization procedure is also the 
simplest when using the EDF Eq. (3), we did not con- 
sider explicitly any other forms. Moreover, nobody ever 
suggested other forms for the contribution of the pair- 
ing correlations to the energy density in the past, as far 
as we are aware. We do not expect however any quali- 
tative changes in our results resulting from considering 
other forms for the superfluid contribution to the energy 
density. 

We looked for a self-consistent solution of the 
Bogoliubov-de Gennes equations corresponding to the 
EDF in Rcl. (3), with a vortex in the pairing field along 
the Oz-axis and Ti/2 of net angular momentum per parti- 
cle, namely A(r) = A(p) exp(i^), where r = (p, <fi, z) are 
cylindrical coordinates. The mean-field depends only on 



2 




FIG. 1. The density profile of a vortex with the symme- 
try axis Oz as a function of kFp, where kF is the value of 
the Fermi wave vector in the homogeneous phase rioo is the 
asymptotic value of the density. The solid line is for the pa- 
rameter set I, while the dashed line is for the set II. 
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FIG. 2. The pairing field profile, where Aoo is the asymp- 
totic value of the pairing field and the rest of the notations 
are identical to those in Fig. 1. 

the radial coordinate p. The quasiparticle wave functions 
have the structure 



/ u a (r) \ _ / u s (p) exp [in<j> + ikz] 

V v a (r) ) ~ \ v s (p) exp [i(n - l)4> + ikz] 




where s labels the radial part of the quasiparticle states, 
n is an integer, k is the wave vector of the quasi- 
particle state along the vortex symmetry axis Oz and 
a = (s,k,n). The presence of the vortex implies a net 
(super) flow of the Fermion superfluid around the vortex 
axis with the velocity profile 

V s (p)e = — i-^ w ;(r)ie J-w Q (r), (5) 

mp n(pj oq> 

where = (y,-x,0)/p, 

The most salient feature emerging from this analysis 
is the behavior of the density as a function of the ra- 
dial coordinate p, the prominent density depletion on 
the vortex axis, see Fig. 1. The density at the vortex 



FIG. 3. The velocity profile in units of the Fermi velocity 
in the homogeneous phase vf = TikF /m and the rest of the 
notations are identical to those in Fig. 1. The dotted line cor- 
responds to the ideal vortex velocity profile V v (p) = h/2mp. 

core is lowered to a value of s=s 0.3poc, while the pair- 
ing field vanishes at the vortex axis as expected, see Fig. 
2. By comparing the actual flow profile with that of an 
ideal vortex, see Fig. 3, one can see that only the frac- 
tion f s (p) = V s (p)/V v (p) of the system is superfluid at a 
distance p from the vortex axis. We have performed also 
calculations by imposing the artificial constraint that the 
pairing field in the homogeneous phase has a larger value 
than the one expected from the self-consistent solution, 
corresponding to the EDF Eq. (3). In that case the den- 
sity depletion at the vortex core becames even larger. In 
hindsight this result could have been expected. Large 
values of the pairing field correspond to the formation of 
atom pairs/dimers of relatively small sizes. When these 
dimers are relatively strongly bound and when they are 
also widely separated from one another, they undergo 
a Bose-Einstein condensation. For a vortex state in a 
100% BEC system the density at the vortex axis van- 
ishes identically. Therefore, by increasing the strength 
of the two-particle interaction, the Fermion system sim- 
ply approaches more and more an ideal BEC system, for 
which a density depletion of the vortex core is expected. 

The value of the coupling constant, and correspond- 
ingly the value of the scattering length can nowadays be 
routinely controlled in experiments [21], following earlier 
theoretical suggestions [22]. One can then expect that it 
would be feasible to create a BCS state as in the exper- 
iment of Ref. [6], stir the system so as to create one or 
several vortices in a manner similar to the experiments 
with dilute Bose gases [1,2,4] and change the ambient 
magnetic field so that the system is brought close to the 
Feshbach resonance. The first quantized vortex will ap- 
pear only at some critical angular velocity, if the system 
is gradually spun up. The vortex quantization could be 
demonstrated for example as in Ref. [1] . At zero temper- 
ature, the single vortex state is the lowest energy state 
with the total angular momentum Nh/2 (here N is the 
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total number of particles) and its decay into lower en- 
ergy states with total angular momenta differing by a 
few number of h/2 is strongly suppressed. The next vor- 
tex will appear at a larger critical angular velocity, and 
so forth. The density profiles of the vortices now develop 
significant density depletions at the cores, which could 
hopefully be subsequently visualized, if they survive the 
expansion when the trap is removed. As discussed in Ref. 
[7] , a fluid with a power law equation of state (the present 
case) will have a simple scaling expansion, see also Refs. 
[5,23]. As a consequence, the density depletion along the 
vortex core should survive the atomic cloud expansion 
upon trap release. Note that a system in which the scat- 
tering length is infinite is highly unusual, as, even when it 
expands and its density decreases, the relative role of the 
interaction does not change, unlike any other interacting 
system, and such a system remains strongly coupled at 
all densities. 

By increasing further the strength of the interaction, 
two atoms will form a bound state. In the limit when 
the size of this bound state becomes significantly smaller 
then the mean interparticle separation the system be- 
comes a BEC of atom pairs. These pairs repel each 
other with an estimated scattering length sa 0.6 ... 2a 
[12,16]. The BEC phase can in principle be described as 
well within the framework of the same formalism outline 
above. Such an approach would be unreasonable, since 
an accurate description of the relatively strongly bound 
two-atom system in terms of extended mean-field quasi- 
particle states would require an exceedingly large num- 
ber of quasiparticle states. The ratio of the quantum 
states required for a reasonable accuracy to the num- 
ber of particles per particle would be at least of order 
l/(na 3 ) 3> 1. The natural description in this limit is in 
terms of bosonic dimers, which weakly repel each other. 
The magnitude of the coherence length in the BEC phase 
is exceeding considerably the mean separation between 
pairs, i=a ^/na ^> 1/n 1 / 3 . This coherence length £b de- 
termines the radius of the vortex core. Thus along with 
the two-atom scattering length a, the size of the vortex 
core could be controlled as well by means of the Feshbach 
resonance. 

In summary, we have shown that the vortex in a dilute 
atomic supcrfluid Fermi gas interacting with an infinite 
scattering length develops a significant density depletion 
along its core, which could be visualized after expanding 
the atomic cloud. 
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